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Abstract 

We present a new method to study disordered systems in the low temperature 
limit. The method uses the replicated Hamiltonian. It studies the saddle points 
of this Hamiltonian and shows how the various saddle point contributions can be 
resummed in order to obtain the scaling behaviour at low temperatures. In a large 
class of strongly disordered systems, it is necessary to include saddle points of the 
Hamiltonian which break the replica symmetry in a vector sector, as opposed to 
the usual matrix sector breaking of spin glass mean field theory. 

LPTENS preprint 96-62 

1 Introduction 

The use of the replica method has turned out to be very efficient in some disordered 
systems. It allows for a detailed characterization of the low temperature phase at least at 
the mean field level. In all the mean field spin glass like problems where one can expect 
the mean field theory to be exact, the Parisi scheme of replica symmetry breaking [JTJ] is 
successful, and at the moment there is no counterexample showing that it does not work. 
On the other hand, the low temperature phase of these systems is complicated enough, 
even at the mean field level. One might hope that the very low temperature limit could 
be easier to analyse, while its physical content should be basically the same. This very 
low temperature limit is also an extreme case where one might hope to get a better 
understanding of the finite dimensional problem. At first sight the low temperature 
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limit is indeed simpler since the partition function could be analysed at the level of a 
saddle point approximation. However it is easy to see that generically this limit does 
not commute with the limit of the number of replicas going to zero. There is a very 
basic origin to this non commutation, namely the fact that there still exist, even at 
zero temperature, sample to sample fluctuations. In this paper we try to develop a 
method of summation over all saddle point in replica space, in order to get the low 
temperature behaviour of glassy systems. The main aim of this paper is to propose 
this new method. We have tested it on some elementary problems which can be solved 
directly. As for its application to more difficult problems, we have also obtained some 
very good approximation to the zero temperature fluctuations a particle in a random 
medium, as well as some interesting scaling relations in the random field Ising model. 
Sect. 2 presents the method and illustrates it on a variety of zero dimensional problems. 
In section 3 we discuss the case of directed polymers in random media with long range 
interactions, where we rederive the scaling exponents using this new method. In sect. 4 we 
study the D dimensional random field Ising model. Perspectives are briefly summarized 
in sect. 5. 



To demonstrate in the simplest terms how the proposed procedure works, we consider 
first some trivial zero dimensional problems. The simplest example is one Ising spin 
a = ±1 in a random field h. The Hamiltonian is: 



2 Zero-dimensional systems 



2.1 The Ising Model 



H = ah 



(2.1) 



where the distribution for the random field is Gaussian: 




(2.2) 



The free energy is: 




(2.3) 



where Dx is the centered gaussian measure of width 
particular, in the zero temperature limit one finds: 
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Let us consider now how this "problem" can be solved in terms of the replica ap- 
proach: 



-PF{ho; (3) = lim 



n->0 



lim 



n->0 



E W}=± iex P {|/3 2 ^(EL l( Ta) 2 }-l 



(2.5) 
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n—>0~, 



n 



■exp{lf3 2 h 2 (2k 
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'k=0 k\(n-k)\ 

In view of the application of the method to more complicated problems we want to 
compute the behaviour at low temperature. This cannot be done naively from a saddle 
point evaluation of the sum at large f3, because of the non commutativity of the limits 
(3 — > oo and n — > 0. Instead we proceed as follows. The term k = 0, which is the 
contribution from the 'replica symmetric (RS) configuration' a a = +1, is singled out; its 
contribution is equal to l + 0(n 2 ), which cancels the (—1) in eq.( [2.5| ). The contributions 
of the rest of the terms (which could be interpreted as corresponding to the states with 
"replica symmetry breaking" (RSB) in the replica vector {cx a }) can be represented as 
follows: 



F(h-B)- lim 1 V ----- r 1 , 

F(h Q , P) - - Inn - ^ T{k + mn 



j—expi-phZQk-ny} (2.6) 
n to oo since the gamma 



Here the summation over k can be extended beyond k 
function is equal to infinity at negative integers. 

Now we perform the analytic continuation n — > 0, Using the relation: 



I> + 1} 
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r(fc + i)r(n - k + v 

Thus, for the free energy (p.6|) one obtains: 



|rwO 



n- 



(2.7) 



-PF(ho; (3) = J2T=i iz ¥ Zi exp(/3 2 ^F) = /+» Dx \ n [i + e xp(2/3/ io x)] (2.8) 



We see that this result coincides with the one fl2.3|) obtained by the direct calculation. 



This is of course no surprise since we have just done an exact replica computation. 
But it exemplifies some of the steps that we shall need below, in particular the proper 
definition and computation of the divergent series appearing in ( |2.8| ) through an integral 
representation. 

2.2 The "Soft" Ising Model 

Consider now the "soft" version of the Ising model described by the double- well Hamil- 
tonian: 



3 



H=--r<f + -^-H (2.9) 

where the random field is described by the Gaussian distribution ( |2.2[ ). We concentrate 
again on the zero temperature limit. Besides, we assume that the typical value of the 
field ho is small (h <C r 3 / 2 ). In this case the field will not destroy the double- well 
structure of the Hamiltonian ( |2.9| ), and (at T — > 0) the system must be equivalent to 
the discrete Ising model considered before. (The "opposite limit" of the random field 
Hamiltonian with only one ground state will be considered in Sec. 2. 3). 

The direct calculation of the zero-temperature free energy is trivial. For a given value 
of the field h <C r 3//2 the ground states of the Hamiltonian ( |2.9| ) are: fa ~ +v^~ + h/2r, 
for h > 0; and fa ~ — v^~ + h/2r, for h < 0. In both cases the corresponding energy is 
Eg{h) ~ — |r 2 — \h\^/r. Thus, the zero-temperature averaged free energy is: 



, 1 9 ,- f +0 ° dh , . h 2 . 1 9 Ih^Jr 
F(h Q ) ~ —r 2 - 2^ / -—fcoq, ) = --r 2 - -^4- (2.10) 



Consider now how this result can be obtained in terms of replicas. The replica 
Hamiltonian and the corresponding saddle-point equations are: 

in i n -r n 

Hn = —rY,^ + - A T,^- >o(£ 0a) 2 (2.11) 

Z a=l 4 a=l Z a=l 

n 

-T<f> a + <f>l = PhlC£<t>a) (2.12) 
a=l 

The "replica-symmetric" solution of these equations (in the limit n — ■> 0) is: a = 
0rs — V^~- The corresponding energy is E rs = — ^nr 2 . This solution (in the limit n — > 0) 
does not involve the contribution from the random field. 

Proceeding along the lines of the Section 2.1, we have to look also for the solution 
of the eqs.( p.l2j ) which would involve the "replica symmetry breaking" in the replica 
vector {fa}: 

a ={ ^ a= }>";>* (2.13) 
r [02 for a = k + 1, ...,n v ; 

In terms of this Ansatz in the limit n — > the replica summations can be performed 
according to the following simple rule: Z)"=i 0a, = kfa + (n — k)(f> 2 — > k(fa — fa). The 
saddle-point eqs.( 2.12 ) then turn into two equations for fa and fa: 

-Tfa, 2 + <Pl 2 = Pkh 2 (fa-fa) (2.14) 
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Assuming that f3kti^ <C r (the explanation of this strange assumption -considering that 
we are interested in the f3 —>■ oo limit!- will be given below), in the leading order one 
gets: 

0!~+Vr; 02 ^ ~Vr (2.15) 



From the eq. (|2.11| ) one obtains the corresponding energy of the above "RSB" saddle- 
point solution: 

E k = -\kT{(j>l - 2 ) + \k{cj>{ - 0!) - \(5h 2 k 2 {^ - 2 ) 2 ~ 

(2.16) 

~ -2(3k 2 h 2 r + O(h 4 ) 

Now, similarly to the calculations of Sec. 2.1 for the zero-temperature free energy one 
obtains, summing the contributions from all these saddle points: 

F(h ) = -lim^o^n-l) ~ 

— hm n ^o ~kl{Z' , rs" ~ 1) — lim ra ^ ^^"rsb" — 

(2.17) 

-hm_ ^-[exp(i^r 2 ) - 1] - I 1 exp(2(3 2 k 2 h 2 r) = 

-P - \ /-« -jt exp(-|x 2 ) ln[l + exp(20 W?)] 



4 ' /3 J-oo ^ ^ 2 

Taking the limit (3 — > oo one finally gets the result 



. . > 1 9 1 /- /" +0 ° cb , 1 9 , 1 9 2/i 0v /r 

which coincides with eq.( |2.10| ). 

It is worth to note that the summation of the series in eq. (|2.17| ) can also be performed 
in the other way: 



F rsb = -\ f) exp(2/3 2 fc 2 / i 2 r) = -L / — ^— exp^V/^r) (2.19) 

K 2«p Z Sin(7T2) 

where the integration goes over the contour in the complex plane shown in Fig. la. 
Then we can move the contour to the position shown in Fig. lb, and after the change of 
integration variable: 

z -> [2f3 2 hlrY l / 2 ix (2.20) 
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in the limit (3 — > oo we have: 

sin(Trz) ~ ^-m[2hlT}- l/2 x (2.21) 

Then, taking into account also the contribution from the pole at x = for the integral 
in eq. (|2.19|) we get: 



rsb 



h y/2r f+°° dx 2 2h 0y /r 



[exp(-x 2 ) - 1] = (2.22) 



27T J-oo X 2 V2VT 



which again coincides with the results ( |2.10|) and ( |2 . 1 8 ) . 



This little exercise with the integral representation of the divergent series in eq. (|2.19|) 
shows in particular that the "effective" value of the parameter /3k — > (3z which enter 
into the saddle-point equations (|2.14|) scales (according to (|2.20|) ) as {hQ^/r)^ 1 . That is 



why in the zero temperature limit the "effective" value of the factor fikti^ ~ h /y/r in 
the eq.( [2.14| ) can be assumed to be small compared to r (for small fields h <C r 3 / 2 ). 

Replica fluctuations 

Because of the non commutativity of the limits n — > and (3 — > oo, one cannot get the 
exact result by keeping only the saddle-point states of the replica Hamiltonian. Actually, 
averaging over quenched disorder involves the effects of sample to sample fluctuations 
which in terms of the replica formalism manifest themselves as the contribution from 
the replica fluctuations. In other words, to get exact result in terms of replicas the 
contribution from the saddle-points is not enough, and one has to integrate over replica 
fluctuations even in the zero-temperature limit. 

This phenomenon can be easily demonstrated for the above example of the "soft" 
Ising model. Les us take into account the contribution from the Gaussian replica fluc- 
tuations near the "replica-symmetric" saddle-point (p a = <fi rs = y/r: 



'rs 



+ (2.23) 
From the eq. (|2.11|) for the "replica-symmetric" part of the partition function we get: 

Z» rs " = exp(|/3nr 2 ) / al^ a exp{-(3 Ea,b( T $ab - \(3hl)^ a ^ b } 

*< exp{i/3nr 2 + - |nln(^r)} ^ 

Therefore, in the zero-temperature limit one obtains the following contribution to the 
free energy: 

1 h 2 

F» Pa » = --r 2 - ^ (2.25) 
4 4t V ' 
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We see that at T = there exists a finite contribution ~ h^/r due to the replica 
fluctuations. In the particular example considered the value of h Q was assumed to be 
small, and this contribution can be treated as a small correction. However, we should 
keep in mind that the contribution from the replica fluctuations in general could appear 
to be of the same order as that from the saddle points. Therefore, the calculations 
we are going to perform in next sections for less trivial examples taking into account 
only saddle-point states can not pretend to give exact results, giving only the scaling 
dependence from the parameters of a model. 



Saddle points 



In the above calculations of the free energy for the "soft" Ising system we have taken 
into account only the contribution from the two minima of the double- well potential. The 
existence of the third saddle-point, which is the maximum at cf) = 0, has been ignored. 
In this particular example such an algorithm looks natural. However, in less trivial 
systems very often it is not easy to distinguish the types of the saddle points involved. 
Moreover, it could be very hard to impose a simple and robust "discrimination" rule 
with respect to different types of saddle-points, which would not block the calculations 
at the very start. 

Because of this, we would like to propose a somewhat modified scheme of calculations 
which takes into account all saddle points. In the above example of the "soft" Ising model 
the third saddle-point (the maximum) is at = 0. Then, instead of the Ansatz ( 2.13Q 
let us represent the replica vector <p a as follows: 

+y/r for a = 1, k 

-y/r for a = k + l,...,k + l (2.26) 
for a = k + I + 1, n 



For the corresponding "energy" (in the limit n — > 0) from the replica Hamiltonian ( 2.11| ) 
one easily finds: 

H kl = ~r 2 {k + 1) - ~ph 2 r(k - I) 2 + O(h 4 ) (2.27) 

Note that in terms of the Ansatz ( |2.26j ) the "replica symmetric" state (k — I — 0), 
(p a = 0o = has zero energy, so that it gives no contribution to the free energy. 
The combinatoric factor in the n — > limit is now: 

„, B (-l)^i( t + pi 



k\ l\(n — k — l)\ k + l kill 

Thus, for the free energy (for (3 — *> oo) we obtain: 
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F(ho) 



This series can be summed up in a similar way as the ones in eqs.( |2.8| ) and ( 2.17|) 



F(h c 



1 f +oc ~ 3 3 

— I Dx\n 1 + exp{-r 2 + fiho^/rx} + exp{ — r 2 - ph s/rx} 

p J —00 4 4 



(2.30) 



In the limit 3 — > 00 one finds: 



1 /-+00 
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+ /^ o 0->| 



2/ipy^ 
v 7 ^ 



(2.31) 



Again, we get the correct result. While it might seem at first sight somewhat "magic", at 
least some aspects of this computation can be understood. In the example considered (as 
well as in the further examples to be studied below) the relevant states, which contribute 
to the free energy, have negative energy —E(h). Then, in the low temperature limit the 
partition function of a given sample is Z ~ exp(+3E(h)). Therefore, in the limit 3 — > 00 
the free energy can be represented with exponential accuracy as follows: 



F(h c 



>(1 + Z) 



1 v^oo 



(-1)" 



(2.32) 



One can easily check that after averaging Z m = Z m and taking into account the contri- 
butions of the two minima of the corresponding replica Hamiltonian H m one recovers 
the series in eq.( [2.29| ). 

The only "magic" rule which should be followed in the direct replica calculations 
is that the "background" state, 0o = (the one with zero energy) in the Ansatz for 
the replica vector 4> a of the type ( |2.26|) should be placed in the last group of replicas. 
Using this rule, the series obtained for the free energy will correspond to the above 



interpretation ( |2.32| ) 



2.3 One- Well Potential 

Consider now how the method works in the case where the Hamiltonian has only one 
minimum: 



H = -4> a - h<p 
a 



(2.33) 
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where <fi > and a > 2, and the random field h is again described by the Gaussian 
distribution ( |2.2j ). For a = 4 this system can be interpreted as the variant of the 
Hamiltonian ( |2.9| ) in the limit of strong magnetic fields. 

In the zero-temperature limit the free energy is defined by the ground state 4>(h) = 
h i/{a-i) for /i > o, and = for /i < 0. Its energy is E(h) = ~^-h a ^ a ~^ for h > 0, 
and E = for h < 0. Therefore, for the averaged zero-temperature free energy we find: 



^ ) = -^ir^h^eM-& o } 



(2.34) 



-^o" 1 x ^ Jo + °° i/k x ~ expHa; 2 ) = A(a) x 
In terms of replicas, the corresponding replicated Hamiltonian is: 

-in i n 

H n =-Y,<K-M £ ^ (2.35) 

a a=l 1 a,b=l 

This Hamiltonian has a trivial "background" extremum at <fi = with zero energy. 
Therefore, following the scheme proposed in the previous subsection, we look for non- 
trivial saddle-point solutions in terms of the following Ansatz: 

<Pa=it ^ a= }>";>* (2-36) 
r [0 for a = k + 1, ...,n v ; 

For the corresponding Hamiltonian and the saddle-point equation (in the limit n —>■ 0) 
one gets: 

H k = -k<f> a - -/tajftV (2.37) 
a I 

a-1 - (3h 2 k(f) = (2.38) 
The solution on this equation and the corresponding energy are: 

<f> = (J3khl)£z (2.39) 

(Note, that although one can try with more " RSB" steps in the replica vector <p a it can 
be easily proved that there exists only one type of the non-trivial solution given by the 
Ansatz ( |2.36| ) ). Then, in terms of the procedure described above for the free energy we 
have: 

1 ~ C — l") fc — 1 a - 2 2(a-l) 

F(h) = -p^-j—eM^m^K- 2 } (2-41) 
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The summation of this series can be performed in terms of the integral representation 
eq.CT): 



I r dz a — 2 2(0.-11 2a 

F(h ) = — R / — ^exp{^— - (Pz^hS- 2 } (2-42) 
2ip Jc zsm{Trz) 2a 

where the integration goes over the contour in the complex plane shown in Fig. la. Then, 
again, we move the contour to the position shown in Fig. lb and redefine the integration 
variable: 

1 2L_ 

z^-h a -Hx (2.43) 

In the limit f3 — > oo we have: 

1 «_ 

sin(7T2;) ~ -rnr/t Q_1 x (2.44) 

and 

a 1 r Ht cu — 9 2(a-i) a 

F(h ) = -hS- 1 - ^ exp{^— (zx)^ 1 } = S(a) x ^ (2.45) 

Thus, we have obtained the correct scaling of the free energy as a function of ho. Note 
however, that although it is also possible to calculate the value of the coeflicien B(a) 
in the integral ( |2.45| ), such a calculation would not make much sense because to obtain 
the correct coefficient (which is given by the integral in ( |2.34| )) one would need to take 
into account replica fluctuations which we have neglected here. 



2.4 The Toy Model 

Let us consider now a slightly less trivial example of a zero-dimensional system which 
cannot be solved by elementary algebra. This system, generally called the "toy model", 
consists of a single degree of freedom evolving in an energy landscape which is the 
sum of a quadratic well and a Brownian potential. The Hamiltonian is: 

H = + V{cj>) (2.46) 
where V(<p) is the random potential described by the Gaussian distribution: 

P[Vm ~ exp{-l / d<j>& 2 } (2.47) 
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The V distribution is characterized by its first two moments: 



- V(0')) 2 = U<t>-<t>' 



V((j>) = 0; (2.48) 



where C is an irrelevant constant This problem was introduced originally as a toy, 
zero dimensional version of the interface in the random field Ising model ||. It has 
the virtue of showing explicitely how the most standard field theoretic methods like 
perturbation theory, iteration methods or supersymmetry get fooled in this problem, in 
the low /i, low temperature limit, by the existence of many metastable states [Q, ||, ^ ||. 
The main point is that at low enough temperatures the usual perturbation theory does 
not work and a qualitatively reasonable theory must involve the effects of the replica 
symmetry breaking. This has been demonstrated within the replica gaussian variational 
approximation [|, |T0|j . 

One quantity which one would like to calculate in such a system is the value of (0 2 ) 
in the limit of the zero temperature. Using simple energy arguments one can easily 
estimate what must be the scaling dependence of this quantity on the parameters /i 
and g. For a given value of <ft the loss of energy due to the attractive potential in 
the Hamiltonian (|2.46|) is ~ fi<ft 2 . Possible gain of energy due to the random potential 
according to statistics (|2.47|) - (|2~4"8"D can be estimated as ~ sfg~\f$- Optimizing the total 
energy E ~ /i</> 2 — ^Jg\/4> with respect to <fi one finds that 

2/3 

W) = C 2 fa (2.49) 

This result tells that the typical energy minimum of the Hamiltonian ( |2.46| ) lies at a 
finite distance from the origin. The scaling ( |2.49| ), which is obviously right, is not so 



easy to derive from some field theoretic methods which could be also used in higher 
dimension, and there is no known exact result for the constant C*2 at the moment. 

Let us try to calculate the value of (0 2 ) in the zero temperature limit using the 
method considered above. The replicated Hamiltonian of the system (|2.46|) is: 



in i n 

H n = -vY,4>l+ -Jg £ \cf> a - b | (2.50) 



2 ^ a 2 

Z a=l z a,b=l 



The corresponding saddle-point equations are: 

n 

n<t> a + (3gY, Si 9 n (<t ) «- ( t ) b) = ® ( 2 - 51 ) 



6=1 
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(Note that in this formula, whenever there is some ambiguity, one should always 
assume that there is at some intermediate step a short scale regularization. Therefore, 
one must interpret for instance Sign(0) = 0.) Let us first look for non-trivial solutions 
of the eqs. ( 12.51]) . It can be easily proven that within the "one-step" RSB Ansatz ( |2.13| ) 
there exist no non-trivial solutions. Let us consider the "two-steps" Ansatz for the 
replica vector a : 



pi for a = 1 , . . . , k 
p2 for a = k + 1, k + I 
/> 3 for a = k + I + 1, ...,n 

From eqs.( |2.5lD one finds the following equations for 01,2,3 (in the limit n 



(2.52) 



0): 



^2) - f3g(k + l)Sign(4>i - 3 ) 
(f>i) ~ (3g(k + l)Sign((f)2 - fa) 
01 ) + f3glSign((j) 3 - 2 ) = 



/x0i + (3glSign{4>i - 

H(j>i + /3gkSign(4>2 

Ufa + PgkSign(4> 3 
The solution of these equations is: 

0i = --(3k ; 

and the corresponding energy is (in the limit n —>■ 0): 





= 



(2.53) 



^ 2 = ; 



-pa - k) 



E, 



hi 



2/i 



-kl(k + l) 



(2.54) 



(2.55) 



It can be proven that there exist no other solutions of the saddle-point equations 
( |2.51| ) with a number of RSB steps larger than two. 

Therefore, (after taking the limit n — > 0) for the "RSB" part of the free energy we 
get the following series (see eq.(|2.2 



rsb 



1_ 

fin 



exp(—j3E, 



Mi 



f3n ^k+l=l kW.(n-k-l)\ 

C.\p{A/.7(// • /)} 



_J_Y-oo (-l)fc+'~i (k+l) 



(2.56) 



where 



k+i 



A 



3^2 



2/i 



(2.57) 



We again carry the summation of the asymptotic series ( |2.56j ) with the integral 
method mentioned in section 2.2: 



12 



rsb 



dzidz 2 T(zi + z 2 + l) r , , \ -| 

■expjA^i 22(21 + z 2 )) 



(3(2i) 2 J Jc {z 1 + z 2 ) sin^O sin(7r2 2 ) r(zi + l)Y{z 2 + 1) 

(2.58) 

where the integrations over z\ i2 both go around the contour in the complex plain shown 
in Fig. la. 

Shifting the contour of integration to the position shown in Figl.b, and redefining 
the integration variables: Z\ 2 — > \ ' %X\ 2 in the limit (3 -»• 00 (A~ 1/3 0) one gets: 

IA 1 / 3 f f f +ao , , r sin(xix 2 (xi + x 2 )) sin(a:i3; 2 (xi -x 2 )) 

F rsb = -^TT^i / ozkW 7 j ^ + 7 \—\) 2 - 59 

p 2ir z J Jo XiX 2 (x\ + x 2 ) XiX 2 [X\ — x 2 ) 

Substituting the value of A = j3 3 g 2 /2fi we finally get the result for the zero-temperature 
free energy: 

?/> V3IX1/6) 

" ^73 47T3/2 (2 ' 60) 

To this piece we must now add the replica symmetric contribution. The saddle point 
equations have the trivial solution: <p a = with the corresponding energy E = H n [cf) a = 
0] = 0. As we want to get a quantitative result for the constant C 2 , we must also include 
the contribution from the replica fluctuations around this saddle point. This cannot be 
done just at the level of integrating the quadratic fluctuations. We shall rather make 
the following (strong) assumption, namely that this whole 'RS' part of the free energy, 
including the replica fluctuations, is given by the Gaussian replica variationnal method 
[|, |S], . We do not have a very convincing argument to support this hypothesis; we 
just point out that this gaussian variationnal method involves the Gaussian integration 
over replica fields which in a sense is "symmetric" with respect to the point <p a = 0. In 
the end the hypothesis is best supported by the good result one gets for C 2 . We denote 
the gausian variational contribution by F rv , and our conjecture is that F = F rv + F rs b. 

According to eq.( |2.46| ): 

W) = 2^=W) - Q — r(V6) (261) 



Using the result of [[J for the value of {4> 2 ) rv we finally get: 

o 2 /3 / 3 IY1/6) \ o 2/3 

(02) = 9 _ IV o 1,00181^ (2.62) 

Wl /i 4 /3 ^( 4vr )i/3 2x/3tt 3 / 2 y /W 3 V ; 

We have compared this result with some numerical simulations of the problem. The 
scaling in fi and g is obviously correct, the only point to check is the pref actor C 2 . 
Choosing for instance the values of the parameters /i = 1 and g = 2^/H (when the replica 
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variational method gives (4> 2 ) rv = 3 ) we obtain from (|2.62|) the analytical prediction: 
(<ft 2 ) ~ 2.3291. The numerical simulation was done at zero temperature, with these 
same values of fi and g. The <fi interval [—8, 8] is discretized in 2N points, on which 
one generates a random potential as in ( |2.46| ). The exhaustive scan gives the minimum, 
from which one computes (</> 2 ). We average over 100000 samples. The number of points 
2N ranged from 2 s to 2 16 , in this regime there is no systematic N dependance. There 
is no systematic error due to the finite width of the interval since we have checked that, 
within our statistics, there is no sample for which the minimum is found with \<f>\ > 7. 
The result of the simulation is ((f) 2 ) ~ 2.45 ± .02. The value predicted by our replica 
saddle point summation is rather close, although there is a clear small discrepancy. 



The result for (0 4 ) 



To be sure that this relatively good agreement of our prediction for (<p 2 ) with the 
numerical result is not just a coincidence we have performed similar calculations for the 
next order correlator (0 4 ). The computations, which are similar to the ones we have 
just presented but more cumbersome, are given in the appendix. The result is: 



For the values of the parameters fi = 1 and g = 2^/v: (when the replica variational 
method gives (4> A ) rv = 27 ) we obtain: (0 4 ) ~ 16.25. The numerical result is obtained 
with the same procedure as above and gives (0 4 ) ~ 17.05 ± .2. Again these number are 
close but there is a significative difference. 

Clearly, the vector type of rsb that we have been using on all these zero dimensional 
problems is somewhat strange, and we cannot assert that we control all of its aspects (in 
particular the fact that the replica fluctuations around the rs saddle point are summed 
by the gaussian variational method is still unclear). However in all these cases, and 
in particular in the non-trivial case of the toy model, we have obtained good results 
using this simple receipe. Therefore we now turn to its application to more elaborate 
problems, starting with systems in one dimension. 



3 Directed Polymers in Random Media 

The problem of a directed polymer in a random medium is an important problem which 
has been much studied recently []n}| . Although the situation in 1+1 dimension, with a 
delta correlated potential, is relatively well understood, there are still a lot of uncertain- 
ties about more complicated cases. 
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We shall consider a one dimensional case with long range correlations of the potential. 
It is described by a one-dimensional scalar field system with the following Hamiltonian: 



where the random potentials V(x, <fi) are described by the Gaussian distribution with 
non-local correlations with respect to the fields 0: 



V(x, <f))V{x', <p') = 8{x - x') [const -g{4>- <f)'y a ] (3.2) 

where < a < 2. 

This problem naturally arises, with a = 1/2, when one considers an interface in the 
two dimensional random field Ising model at low temperatures: then the field cf) just 
describes the lateral fluctuations in the interface, in a solid on solid approximation. 

One first basic question that we would like to answer concerns the scaling behaviour 
of the lateral fluctuations. Let the value of the field 4>(x) be sticked to zero at the 
origine: <p(x = 0) = 0. Then one would like to know how the average value of the field 
at x = L, (4>(L) 2 ), scales with L: 



(0(L) 2 ) = U" 1 / dfaffi / D4>(x) exp(-f3H[<f>(x), V]) ~ L ^ (3.3) 

V J J0(O)=O / 

where the partition function Z (for a given realization of the random potential) is given 
by the integration over all the trajectories 4>(x) with only one boundary condition <fi(x = 
0) = 0. The 'wandering exponent' ( has been computed in the case of local correlations 
of the random potential, it is then equal to 2/3 [TjJ. In the case of non local correlations 
such as ( |3.2|) , it is believed that this exponent should be equal to 3/2(2 — a) at small 
enough a. This is the result that is obtained from the gaussian variational Ansatz |JJ], 
and it can also be derived from a mapping to the Burgers (or the KPZ) equation and a 



study of this equation through a dynamical renormalization group procedure [13 



A simple derivation of this scaling can be obtained by an energy balance argument 



a la Imry Ma [14|. Let the value of the field be equal to 0o at x = L. Then the loss 
of the energy due to the gradient term in the Hamiltonian ( |3.1| ) can be estimated as 
Eg ~ (j>o/L. The gain of energy due to the random potential term, according to eq. ( |3.2| ), 
can be estimated as Ey ~ —y/L^fgcf)^. Optimizing E g and Ey with respect to 0o one 
finds: 



3 1 



£2(2-aJ g2(2- a) (3.4) 



In this section we will demonstrate how this result can be obtained in the zero- 
temperature limit in terms of the proposed replica saddle point method. The replicated 
Hamiltonian is: 
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H„ 



dx 



1 " , dcj) a (x) 
dx 



) 2 + ^9 E - 6 (X)) 

Z a,6=l 



2a 



(3.5) 



Strictly speaking, the systematic way of solving this problem following our general 
method is the following: one must find n saddle point trajectories (j) a (x) for fixed n 
boundary conditions 4> a (L), then one has to derive the corresponding energy H n [4> a (L)}, 
and finally one has to find the saddle point solutions with respect to the values of 4> a {L). 

Here we shall follow a much simpler strategy. Since it is obvious that there always 
exists the trivial solution <p{x) = 0, we will suppose that the correct scaling can be 
obtained simply by taking into account one non-trivial saddle-point trajectory. In other 
words, from the very begining we are going to look for the saddle point solutions within 
the following Ansatz: 




for a 
for a 



l,...,k 
k + 1, n 



(3.6) 



Comparing this Ansatz to the zero dimensional exercises of the previous section, we see 
that it should amount to assuming that the lowest lying configuration dominates. This 
is certainly true since one knows |L5] , |l6f that the metastable states have an excitation 
energy which scales as L", with u — 2£ — 1. Substituting this Ansatz into the replica 
Hamiltonian (|3.5| ) in the limit n — > one gets: 



k dx 
Jo 



1 d(j)(x) 2 
2 { dx ' 



f3kg(f) 



2a , 



X) 



(3.7) 



As usual (see the previous section) the free energy is defined by the series: 



F(L) 



-1 oo 

■-T 



k=l 



k 



exp(-(3H k 



(3.8) 



where the value is defined by the corresponding saddle-point solution for (j)(x). 
The saddle point trajectory is defined by the following differential equation: 



dx 2 



—2a(3kg(j) 



2a- 1 



with the boundary conditions: 0(0) = and 4>(L) 
integrated: 



(3.9) 

b . This equation can be easily 



<t>( x ) d(j) 



xJ2(3kg 



(3.10) 



where the integration constant A is defined by the boundary condition: 
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k - <P' 2a v ' " 

Substituting this solution into the Hamiltonian (|3.7|) , we obtain after some simple alge- 
bra: 



Hk — k 



-(3kg\L + ^2j3kg (*° d<j)y/\ - <p 

J 



(3.12) 



Taking derivative of Hk with respect to 0o (and taking into account the constrain ( |3.11|) ) 
one finds the following saddle-point solution: 

(j) ~ Lii ((3kg)^o^) (3.13) 
and A = 0q q . Its energy (|3.12| ) is: 

(cOnSt) , „, , 2-a r l+a 1 . . 

iZ* = -^-p^W** L— g— (3.14) 
Now we proceed as before, introducing an integral representation of the series ( |3.8| ) and 

i 1+2 1 

a rescaling of the integration variable by a factor -^L 2 ~ a g 2 ~ a . Then we get the scaling 
of the free energy: 

F(L)~L^g^ (3.15) 
from which we obtain the scaling of <po as a function of L: 

(f>o(L) ~ L2(2^)g2(2^) (3.16) 

which coincides with the result ( |3.4|) given by the naive energy arguments, as well as by 
more elaborate calculations. 

Although the example demonstrated in this section provides no new results we hope 
that the proposed method could turn out to be also useful when applied for directed 
polymers with smaller a, or in larger dimension. 



4 Random Field Ising Model in D dimensions 

Since the topic of the random field Ising model covers an enormous amount of litterature 
(see e.g [0), it would be rather difficult to give any brief introductory review. Here, 
however, we are mainly concerned with how the method we have proposed before works 
in various situations. Therefore, we will concentrate only onto one particular aspect of 
the problem. 
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It is well known that the main problem in the studies of the low temperature phase 
in the random field Ising model is that one has to perform the summation over numerous 
local minima states, which seems to be impossible to do within the framework of the 
usual perturbation theory It has been proposed recently that, because of these local 
minima states a special "intermediate" (separating paramagnetic and ferromagnetic 
phase) spin-glass-like thermodynamic state could set in around the critical point, and 
moreover, this state is characterized by a replica symmetry breaking in the corresponding 
correlation functions [ l8| . At low temperature, and when the width of the distribution 
of the random field is not too small, the same phenomenon must be present. It was 
proposed long ago @, and elaborated later on in ||, that the metastable states in 
this regime should be characterized by some "instanton in replica space" . Our method 
provides one more step in the elaboration of this idea. 

We consider the random field Ising model in terms of the usual Ginzburg-Landau 
Hamiltonian in D dimensions: 



H 



d D x 



I(V0) 2 + ^ 2 + ^# 4 -^)^ 



(4.1) 



where the random fields h(x) are described by the 5-correlated Gaussian distribution: 



p[h(x)) = n 



exp 




(4.2) 



The corresponding replica Hamiltonian is: 



H„ 



d D x 



in in in 

Z o=l Z 0=1 ^ a=l 



-, h l £ 0a< 



a,6=l 



(4.3) 



According to the procedure developed in previous sections we are going to look for 
the most simple non-trivial saddle-point solutions at the background of the trivial one, 
<f> a (x) = 0. In terms of the Ansatz: 



[x 




for a 
for a 



l,...,k 



The replica Hamiltonian 

H k = k f d D x 



k + 1, n 
) reads in the n —>■ limit: 



(4.4) 



(4.5) 



Consider for simplicity the situation at r = (Notice that here we work close to the 
critical temperature. The use of our saddle point technique allows to study the system 



at the tree level, which is supposed to give the leading singularities close to T c [|i~9f ). 
The corresponding saddle-point equation is: 
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- A(f) - X(f) + g(j) 3 = (4.6) 
where A = h^k. As usual, the free energy is given by the series: 

oo / -i \k— 1 

F(h ) exp(-tf fc ) (4.7) 

fc=i K 

where the value of Hk is defined by the corresponding saddle-point solution of the 
eq.flU). 

At this stage we see that the situation is getting rather different from the ones 
studied in the previous sections. If we would choose the obvious space-independent 
solution cj) = (A/g) 1 / 2 , we would find that the value of Hk is proportional to the volume 
V of the system: Hk = — \k(\ 2 /g)V = —-^k^hfy. Then, the summation of the series 

( ^4.7| ) would immediately yield a free energy proportional to V 1 / 3 and not to V. Therefore 
this solution, as well as any other solution with an energy Hk proportional to the volume 
of the system, is irrelevant for the bulk properties. 

Thus, we have to look for localized solutions: the ones which are local in space 
(breaking translation invariance) and which have finite energy. Let us assume to start 
with that such an "instanton"-type solution exists (see below), and that for a given k it is 
characterized by the spatial size R(k). Then, if we take into account only one-instanton 
contribution (or in other words we consider a gas of non-interacting instantons), due to 
the trivial entropy factor V/R D (this is the number of positions of the object of the size 
R in the volume V) we get a free energy proportional to the volume: 

k=l 

where must be finite (volume independent). 

It is easy to understand that the equation ( f4.6|) indeed has localized solutions. Let 
us assume that the value of the field <fi{x) is such that A0 2 ^> g<f) 4 . Then in a first 
approximation the saddle-point equation (|4.6|) is linear: 



Acf> + \<p = (4.9) 

The simplest possible spherically-symmetric solutions of this equation in D dimensions 
are the well known Bessel-type functions. In particular there exist oscillating solutions 
which have a finite value 0(r = 0) = 0o at the origin and which decay to zero at r — > oo 
(like ~ r^i - 1 )/ 2 sinr). For example, in dimension D = 3 this solution is simply: 

rv A 

In dimensions D these solutions have a finite spatial scale: 
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R(k) = \- 1/2 = (h 2 k)~^ 2 (4.11) 

and finite energy: 

H k = -(const) fc^A"^ (4.12) 

At the level of the equation (|4.9f) itself, the value of 0o remains arbitrary (the equation is 
linear). On the other hand, from the point of view of the energy this is not an extremum 
since the energy explicitely depends on the value of O (this is the saddle-point solution 
for the fixed boundary condition <p(r = 0) = </>o). If we would let the value of 0o be 
free in the absence of the non-linear term g<f) 4 it would, of course, fall down to infinity. 
However, if we take into account the term gcf) in the "exact" Hamiltonian" ( [4.5| ) it is 
natural to expect that 0o will stabilize around the saddle-point value 

,2 A 



(4.13) 

The above qualitative arguments can be easily verified for the model double-well poten- 



tial: U(<p) = —\4> 2 f° r |0| < y^/g an d U((p) = +oo for \<p\ > yX/g, taken instead of 
the "real" one: U(4>) = —\4> 2 + \g<^- In this case for any |0 O | < y^/g there exists the 



exact Bessel-like saddle-point solution with finite energy ( |4.12j ), and real extremum of 

the Hamiltonian would be achieved at 0o = =t\A/ 9- 

Let us calculate the contribution of such solutions to the free energy. Substituting 
into the series ( |4.8|) the energy of the solution ( (4.12|) , the estimate for the value of 0o 
(|4. 13|) and the characteristic size of the solution ( |4.11| ), together with A = h\k we get: 



\k-i 



OO /_-| . 

F(h ) VJ2 ~ — l (h 2 k)%exp 

k=i K 



(const) 2 6-d 



gh 



2 



(h 2 kY 



(4.14) 



We see that the series is getting strongly divergent only at dimensions D < 6. This is the 
only regime where the considered saddle-point solutions provide a relevant contribution. 

Now, following the scheme developped in the previous sections, we turn to the inte- 
gral representation and rescale the integration variable by a factor (gh'fy'^D fa- 2 ( which 
gives a free energy with the following scaling in the limit gh\ 1: 

F -M- ~ !(,*>) A (4.15) 

Besides, using the same scaling k ~ (gtifys^hQ 2 for the characteristic spatial scale 
of the saddle-point solutions ( f4.11|) , which could be interpreted as a kind of disorder 



induced finite correlation length (near T = T c , as we shall see in more details below), 
we obtain: 
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RM ~ {gK)-*=B (4.16) 
In the same way one gets the estimate for the value of the "disorder parameter" 2 ~ 



l2 *-f„u2 



'0 ~ 



~{9K)^ (4.17) 

Finally, one can easily obtain the estimate for the value of the temperature interval 
t c around T c where all the above qualitative calculations make sense. Formally the 
derivation of the saddle-point solutions has been done for r = 0. Actually, according 
to the replica Hamiltonian ( |4.5| ) the calculations should remain correct until \t\ <C h^k. 
Using the above scale estimate for k one finds the upper bound for the value of r: 



|r|<r c ~(^)— (4.18) 

This value of r c can be interpreted as the estimate for the temperature interval around 
T c where the supposed disorder dominated (spin-glass type) phase sets in. 

Of course, the procedure proposed in this section is still incomplete. In a selfconsis- 
tent approach one should study the effects produced by the interactions between these 
instanton solutions, not talking about the effects of the critical fluctuations. At the 
present stage we are not able to say anything about the ferromagnetic phase transition 
itself and in particular about the behaviour of the corresponding ferromagnetic order 
parameter. 

Nevertheless, we shall now show that these simple replica instanton estimates are 
quite reasonable and can in fact be recovered in terms of (completely independent) simple 
scaling arguments. Indeed, let us come back to the original random field Hamiltonian 
( |4.1|) . Configurations of the field <f)(x) which correspond to local minima satisfy the 
saddle-point equation: 

- A(f)(x) + T<p(x) + g<f> 3 (x) = h(x) (4.19) 

Let us estimate at which spatial and temperature scales the random fields give a dom- 
inant contribution. We consider a large region of linear size L> 1. The spatially 
averaged value of the random field in this region is: 

h(Q L ) = / d D xh(x) (4.20) 

Correspondingly, the typical average value of the random field in this region of size L 
is: 

h L = \W[Ql)} 1/2 = h L- D / 2 (4.21) 
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Then the estimate for the typical value of the order parameter field 4>l in this region 
can be obtained from the saddle-point equation: 



T<pL + 



h, 



Then, as long as: 



(4.22) 



(4.23) 



t<Pl < 

the typical value of 0l inside such clusters is dominated by the random field: 

<Pl ~ {h L /gf" ~ (hy/3 L -D/e (424) 
9 

Now let us estimate up to which characteristic size of the cluster the external fields 
can dominate. According to (|4.23|) and ( f4.24[ ) one gets: 



L « 



-3/D 



(4.25) 



On the other hand, the estimation of the order parameter in terms of the equilibrium 
equation ( |4.22|) can be correct only on length scales much larger than the size of the 
fluctuation region which is equal to the correlation length (of the pure system) R c ~ t~ u . 
Thus, one has the lower bound for L: 



L » t~ v 

Therefore, the region of parameters where the external fields dominate is: 



<< 



M) 1/D 



T 



3/D 



or 



(4.26) 



(4.27) 



_3-i/D 



t" - « gh 2 
Such a region of temperatures near T c exists only if: 

vD < 3 



(4.28) 



(4.29) 



In this case the temperature interval near T c in which the order parameter configurations 
are mainly defined by the random fields is: 



(4.30) 
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In the mean field theory (which correctly describes the phase transition in the pure 
system for D > 4) v — 1/2. Thus, according to the condition fl4.2U| ) the above non- 
trivial temperature interval r c exists only at dimensions D < 6. Substituting v — 1/2 
into ( f4.30| ) we get: 



r c (h ) ~ (gh 2 )— (4.31) 



Then, the random field defined spatial scale can be estimated from (|4.25| ) 



L c (h ) ~ (<7^)~^ (4.32) 

Correspondingly, the typical value of the order parameter field at scales L c (h ) is 
obtained from the eq. ( |4.24j ): 



1 ~ ~{ghl)^ (4.33) 



The energy density is estimated as y ~ 0l c ^l c - Taking into account (|4.21|) and 
(|Q3D we find: 



f ~ -{ghl)^ (4.34) 

We see that we get through these simple arguments a region around T c where the 
disorder induces a finite correlation length. Furthermore the estimates for -p, L c , and 
r c perfectly coincide with the results obtained in terms of our previous replica saddle- 
point method, eqs. (|4.15|) - ([4.18 ). Both approaches clearly hold only in a regime where 



critical fluctuations can be neglected. 



5 Conclusions 

We have proposed a method to analyse random systems by summing up various saddle 
point contributions in the replicated Hamiltonian. We think that it may open a new 
route in this type of study. In particular, the application to finite dimensional systems, 
which we started here with the directed polymer on one hand, and with the random 
field Ising model on the other hand, looks quite interesting. Indeed we have seen on this 
last case how this method allows to take into account instanton contributions which 
are usually out of reach of most analytic methods in these systems. Such instanton 
contributions have been argued to be important for a long time (0, [|). We think we 
can get them under control with the present approach. 

Clearly our method is still not totally understood in all details. We have pointed 
out that it involves one single basic rule, stating the way one has to order the various 
saddle points in replica space. Within this hypothesis it gives reasonable results in all 
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the cases we have checked so far, but of course more studies are needed to justify this 
hypothesis. 



6 Appendix: Computation of the fourth moment in 
the toy model 

Using the saddle point solution Q2.54D we have: 



<0 4 U = E n k+ i=i m n-k-iy Hi + l(f>i + (n - k - 1)4} exp{-(3E kl } - 

(6.1) 

(f ) 4 E&=i ( ~ 1 j£ , " 1( ffi' fc*(fc + 0(3fc 2 + 3/ 2 - 5kl) «p{AW(* + 0} 
Proceeding similarly to the calculations of the free energy F rsfe (|2.58|) - (|2.59|) we get: 

(^)rsb = ) 4 ^{ _ (2^ / ic ( Zl +2 2 )sfn( 1 7S : i)sm(7r22) f(il+l^^+l) ( 3z l + 3 ^2 - 5^i^ 2 ) X 

exp[A2i2 2 (zi + z 2 )]} 

(6.2) 

Shifting contour to the position in Fig. lb and redefining — > A l ^ix\^ in the limit 
/3 -> 00 (A~ 1/3 -> 0) we find: 



(6.3) 

exp[-i\xix 2 (xi + x 2 )}} 

Taking into account the contribution from the pole at = after somewhat painful 
algebra we finally obtain the following result: 

# 4 / 3 17 v / 3[sin(7r/12) + cos(vr/12)] 



^ rsb ^ 3v^T(l/6)sin(7r/6) (6 ' 4) 
Taking into account the contribution from the replica fluctuations [|§ : 

^4/3 27 

= -JK^J^W (6 ' 5) 

for the fourth order correlator we get the final result: 



^4/3 27 ^4/3 i7^/3[ s i n ( 7r /i2) + cos(tt/12)] 



l ) (^rv + i^rsb ^8/3(^)2/3 ^8/3 3^(1/6) Sin(7r/6) (6 ' 6) 
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Figure captions 

Fig.l The contours of integration in the complex plane used for summing the series, 
a) The original contour, b) The deformed contour. 
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